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ABSTRACT 

The extragalactic background light (EBL) from the far infrared through the visible and extending 
into the ultraviolet is thought to be dominated by starlight, either through direct emission or through 
absorption and reradiation by dust. This is the most important energy range for absorbing 7-rays 
from distant sources such as blazars and gamma-ray bursts and producing electron-positron pairs. In 
previous work we presented EBL models in the optical through ultraviolet by consistently taking into 
account the star formation rate (SFR), initial mass function (IMF) and dust extinction, and treating 
stars on the main sequence as blackbodies. This technique is extended to include post-main sequence 
stars and reprocessing of starlight by dust. In our simple model, the total energy absorbed by dust is 
assumed to be re-emitted as three blackbodies in the infrared, one at 40 K representing warm, large 
dust grains, one at 70 K representing hot, small dust grains, and one at 450 K representing polycyclic 
aromatic hydrocarbons. We find our best fit model combining the Hopkins and Beacom SFR using the 
Cole et al. parameterization with the Baldry and Glazebrook IMF agrees with available luminosity 
density data at a variety of redshifts. Our resulting EBL energy density is quite close to the lower 
limits from galaxy counts and in some cases below the lower limits, and agrees fairly well with other 
recent EBL models shortward of about 5 /Ltm. Deabsorbing TeV 7-ray spectra of various blazars with 
our EBL model gives results consistent with simple shock acceleration theory. We also find that the 
universe should be optically thin to 7-rays with energies less than 20 GeV. 

Subject headings: galaxies: active — diffuse radiation — gamma rays: observations — stars: luminosity 
function, mass function — stars: formation 



1. INTRODUCTION 

The extragalactic background light (EBL) from ^ 
10-3 - 10 eV (- 0.1 " 1000 ^m) is thought to be dom- 
inated by starlight, either through direct emission, or 
through absorption and reradiation by dust. This is 
the most important energy range for photons with en- 
ergy El (TeV) interacting with long- wavelength photons 
with and absorbing 7-rays from distant sources, with the 
threshold 77 condition implying that 



Ei{TeV) 



0.26 



{l + z) EEBL{eV) 



The EBL is also an important target radiation field for 
energy- loss and/or photodisintegratio n of the highest en - 
ergy cosmic ray protons and ions fe.g.- lPuget et al.lll976D . 
which results in a source of high-energy neutrinos (e.g. 
IStanev 2 00^. Direct meas urement of the EBL is difficult 
(see lHauser fc Dwek '20011 for a review) due to contami- 
nation by foreground zodiacal and Galactic light. Galaxy 
counts may also be used to estimate the EBL, but the 
unknown number of unresolved sources results in a lower 
limit. The general picture is a component peaking at 
around 1 /im from direct starlight emission and one peak- 
ing at ~ 100 /im from re-emission of absorbed starlight 
by dust. Due to different modeling approaches and un- 
certainties in underlying model parameters, the intensity 
and shape of the EBL spectrum remains controversial. 

A wide range of EBL models have been developed. 
One class of models approaches the problem by using 
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IR data from local galaxies and extrapolating their evo- 
lution to higher redshifts and shorter wavelengths (e.g., 
IMalkan fc Stecke^ [19981 [2OOII: IStecker et al.ll2006D . An- 
other class of models uses semi-analytic merger-tree mod- 
els of galaxy formation to determine the star formation 
history of the universe for a forward-evolution calculation 
(IPrimack et al.l [200l [20051 [20081 IGilmore et all [20081 
I2009D . The class of models considered here focuses on 
the primary energy source of the emission, namely the 
stars. In such models, one integrat es over star forma- 
tion rates and s tellar properties (e.g. Salamon fc S teckeil 
[T99I iDwek et al . 1998; Kneiske et al. 200 2. 2004). Al- 
terna tively, luminosity density data ( Fra nceschini et al.l 
|2008() is used to estimate the EBL energy density and its 
evolution with time. Since the EBL absorbs very high 
energy (VHE; w 0.1 - 100 TeV) 7-rays from blazars, 
one can use this mechanism to put upper limits on 



the EBL energy densit y (e.g. , IStecker fc de Jageil 



Stanev fc Fraiiceschinil [19981: lAharoiiian et al 



1993 



2006 



Mazin fc Rauell2007l:lFinke fc Razzaauell20 09'). Disagree- 



ment about the intrinsic 7-ray s pectra of blazars (e.g., 
IStecker et all l2007t iBottcher et' al. 2008) has resuhed, 
however, in a lack of consensus about the meaning of 
these upper li mits. Detection of a 33 GeV photon from 
GRB 090902B (lAbdo et al.ll2009a|) and a 13 GeV photon 
from GRB 080916C (lAbdo et al.ll2009E) with the Large 
Area Telescope (LAT) onboard the Fermi Gamma Ray 
Space Telescope are used to test different EBL models, 
which m ay be a promising new way to constrain the EBL. 

We ([Razzaque. Dermer. fc Fink3 120091 hereafter 
RDF09) recently developed models for starlight lumi- 
nosity density and EBL intensity by assuming that stars 
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are blackbodies, and integrating over stellar properties 
on the main sequence. Here we extend the model of 
RDF09 to stars that have evolved off the main sequence 
and include re-emission of absorbed starlight by dust, 
noting that the models of RDF09 remain valid for 
the ^ 1 eV range. In § [21 analytic expressions for 
the radiation from stars and re-emission by dust are 
presented. Numerical calculations of the luminosity 
density and EBL energy density are presented in § [3l 
We explore the effects of our EBL model on absorption 
of distant 7-rays (§ U) and conclude with a discussion 
on our results and future research (§[5]). 

2. FORMALISM 

We briefly describe the RDF09 model for background 
starlight and our recent improvements which include the 
emission from post-main sequence stars and dust. In- 
tegrating over star formation in this manner is similar 
to several other models (e.g. iSalamon fc Steckeii [T998I : 
IDwek et al.l[l998l : IKneiske et al.ll2002ll2004D . 

2.1. Direct Starlight Emission 

Stars with dimensionless mass m — M/Mq and age 
are assumed to emit as blackbodies. The photon density 
of a blackbody is given by 



n*(e; m, 9) = 



dN 
dedV 



Stt 



exp[e/e] 



(1) 



where e = hv /mf.c^ is the dimensionless photon energy, 
Ac = h/rrieC « 2.42 x 10~^° cm is the Compton wave- 
length, and = kBT/rrieC^ is the dimensionless effective 
temperature. The total number of photons emitted per 
unit energy per unit time from a star of radius i?(m,t*) 
and effective stellar temperature Q(m,ti,) is 



Ni,{e;m,tiJ) 



dN 
dedt 



= TiR{m, t^fc n^e; e(m, t^) . (2) 



To determine the luminosities and radii of the stars, 
L(m,ti,) and R{m^t^), respectively, as well as the time 
stars spend on various portions of the Hertzsprung- 
Russell diag ram, we used t h e stel lar formulae from the 
appendix of lEggleton et al.l (|1989D . These formulae ap- 
proximate the stellar parameters along the main se- 
quence, the Hertzsprung gap, the giant branch, the hor- 
izontal branch, the asymptotic giant branch, and the 
white dwarf phase for stars of solar metalicity. Thus, we 
assume all stars emitting since star formatio n began have 
solar m etalicity. Note that eqn. (A15) of Eg gleton et al.l 
(jl989( ). which describes the luminosity of the base of the 
giant branch, should be 



L 



2.15AP + 0.22M5 



BGB 



1 + lAx 10-2M2 + 5 X 

(lEggleton et al.l [1990^. Also, eqn. (A22) of 
Eggleton et al.l (|1989i) . which describes the time a 
star spends burning Helium, should be 

tMsLo 



tHe 



Lffe(MO-42 + 0.8) 



(C. Tout and P. Eggleton, private communication). We 
modify eqn. (A27) of lEggleton et al. (,19891 which de- 



scribes the luminosity of a white dwarf in our calcula- 
tions, so that it reads 



L ^ 



40 



{t - t 



WD 



0.1) 



1.4 



in order to avoid a singularity when t = two 
(iHurlev et al.l [2000l) In the above corrections to 
lEggleton et al.l (|l989l ). we use their notation, so that M 
is the star's mass in units of Mq, Lbgb and Lq are in 
units of Lq, and t, tne, tMS, and two are in units of 
Myr. 

Once a star's luminosity and radius have been deter- 
mined, its temperature can be found by 



ksTQ ( L {m,t^) 



1/4 



Rr. 



R{m,ti,) 



(3) 



where Tq = 5777 K is the effective solar temperature, 
Lq = 3.846 X 10'^^ erg s"-^ is the solar luminosity, and 
Rq — 6.96 X 10^° cm is the solar radius. 

We will use unprimed symbols to refer to quantities 
measured in the frame comoving with the Hubble flow, 
and primed or double primed quantities to refer to quan- 
tities in the proper frame of a galaxy at some redshift z. 
Observations made from the solar system, which has a 
small peculiar velocity with respect to the Hubble flow, 
are essentially in the comoving frame. The comoving lu- 
minosity density (i.e., the luminosity per unit comoving 
volume, or the emissivity) of the universe as a function 
of comoving energy e at a certain redshift z (in units of, 
e.g., W Mpc""^) can be found from 



ef"''-'{e;z)^mech^ 



dN 



dt de dV 
dm ^(m) 





dzi 


dt^ 


J z 




dzi 



^Pizi) iV*(e;TO,i*(z,zi)) . (4) 



The luminosity density is dependent on the initial mass 
function (IMF), ^(m), the comoving star formation rate 
(SFR) density (i.e., the star formation rate per unit co- 
moving volume), ip{z), and the fraction of photons which 
escape a galaxy, /esc(e) and avoid being absorbed by in- 
terstellar dust. The relationship between cosmic time 
and redshift is given by 



dt^ 
dz 



Hail + z)^n^{i + zy + n^ 



(5) 



in a flat ACDM cosmology, 
parameters Hq — '^'^ '"^ 
Oa = 0.7. 



We assume cosmological 
= 70 km s-i Mpc-^, n„, = 0.3, and 

A = 0.7. 

iDriver et al.l ( 2008[ ) have applied the dust model 
of iPopescu &: fH^ (poM ) to a survey of - 10^ 
near by galaxies from the Millennium Galaxy Cata- 
log (I Allen et al.l 120061) to determine the wavelength- 
dependent escape fraction of photons in the local uni- 
verse. RDF09 have fit this with a series of power-laws 
and we use this to compute /esc(e)- We also assume that 
any photon with mgC^e > 13.6 eV is absorbed by galactic 
and intergalactic H I gas. These UV photons are not re- 
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processed in our model, and we assume their net energy 
makes a small contribution to the total EBL intensity. 

In this work we choose the limits of integration rrimin = 
0.1, rrimax — 100, and Zmax = 6, although our EBL 
intensities and luminosity density results at low z are not 
strongly dependent on the upper limits. Thus, the model 
of the stellar component does not have any adjustable 
parameters once an IMF and SFR have been chosen. 

To test the accurac y of approximat i ng sta rs as black- 
bodies and the simple 'Eggl eton et al.l (|1989l ) stellar for- 
mulae, we computed the spectra of simple stellar popu- 
lations (SSPs) for various ages from 



"k [fim] 



dN 



he 



dtdX 

dm £^{m) Ni,{e;m,ti,) (6) 



usmg a ISalpeteil (UOSl IMF, i.e., ^(m) oc m , from 
mmin = 0.1 < m < rrimax = 100. Eqn. © approx- 
imates the spectrum of a cluster of stars that were all 
born at exactly the same time, and all now have an age 
of t*. SSPs at high spectral resolution were calculated by 
iSruzual fc CharlotI (|2003l) using various stellar spectral 
libr aries calculated with detail ed stellar structure codes. 
The lBruzual fc CharlotI pQOl SSPs with solar metahc- 
ity are compared with our approximations in Fig. [TJ Our 
approximations agree fairly well longward of ^ 10'^ A (0.1 
/im), particularly for stars between 1 and 10 Gyrs of age. 
This age range of stars when weighted by the SFR dom- 
inates the radiative output between 0.1 and 10 ^m, and 
this is where our approximations are best. Shortward of 
this wavelength the agreement is poor, as large as several 
orders of magnitude. However, at < 0.1 /im most pho- 
tons will be absorbed by H I gas at low z, and at high z 
where stellar populations are younger agreement is much 
better at these wavelengths, differing by no more than 
a factor of ~2-3. We have computed a test EBL model 
using a Salpeter IMF f or both the blackbody approxima- 
tion and using the Bruzual fc C hariot (2003) SSPs and 
found that the difference is no more than « 25%. As 
we show below, this is larger than the difference without 
including post-main sequence stars f? l3.3p . We therefore 
consider our blackbody approximation quite good for cal- 
culating luminosity densities and EBL energy densities, 
particularly compared to other uncertainties in these cal- 
culations. 

2.2. Dust Emission 

The starlight that is absorbed by dust is reradiated 
in the infrared. There are generally thought to be 
three major dust compo nents in the interstellar medium 
(e.g. iDesert et al.lll990[ ): (1) a large grain component 
that absorbs optical light and reradiates in the far-IR, 
found in and around star-forming regions; (2) A small 
grain component that absorbs the far UV and reradi- 
ates in the near-IR, located throughout the disk of spi- 
ral galaxies and is responsible for most of the observed 
dust radiation and (3) a component from polycyclic aro- 
matic hydrocarbons (PAHs) which emit as broad emis- 
sion lines and are not generally i n thermodynamic equi- 
librium with their environment (|Dwek et al.l 119971) . al- 
though they are sometimes modeled as very hot black- 
bodies fe.g. iKneiske et al.,,2002. ). 




A, [A] 

Fig. 1. — SSP spectra for various population ages, for a Salpeter 
IMF, normalized to 1 Mq. The thin lines are the models from 
[Bruzual & Chariot (2003) and the thick lines are our calculations 
assuming stars arc blackbodies. Photons with wavelengths short- 
ward of 912 A (13.6 eV) contribute a negligible amount to the total 
energy budget and are assumed to be completely absorbed. 



We assume the dust emits as a combination of three 
blackbodies, representing three types of dust: a 40 K 
blackbody representing warm, large dust grains, a 70 K 
blackbody representing hot, small dust grains, and a 450 
K blackbody representing P AHs. These t emperatures 
are similar to those used bv iKneiske et al.l (j2002). By 
setting the luminosity density from dust emission equal 
to the luminosity density from starlight absorbed from 
dust, we calculate the infrared emission self-consistently. 
The emissivity absorbed by dust is 



1 



/esc(e) 



[I- fesc{e)]f ie;z) 



where j«*'"'''(£; z) is the stellar emissivity 12. ip . A frac- 
tion of this will be reradiated in each component, so that 



fn / de 



1 



/esc(e) 



[1~ fescie)] f {e;z) 



(7) 



where /„ is the fraction of the absorbed emissivity rera- 
diated in a particular dust component, 0„ = k^Tn / meC^ 
is the dimensionless temperature of the dust component, 
and the subscripts n —\, 2, and 3 refer to the warm dust, 
hot dust, and PAHs, respectively. Since the dust radiates 
as a blackbody, i.e.. 



j«(e;6«) = j° 



exp[e/e„] - 1 



(8) 



the integral in the right side of eqn. (O can readily be 
performed, so that 



dejn(e;e„) = — j°e;| 



(9) 



This allows one to determine the blackbody normaliza- 
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tion, 



15- 



fnjde 



(10) 



Hence, the comoving luminosity density due to dust is 



■dust I 



15 

3 



j de 


1 

Jesc(e) 


fn 


£4 



64 exp[e/e„] - 1 



(11) 



We assume that fescj fm and 0„ do not vary with z. The 
values for /„ and 0„ were chosen by trial and error to fit 
the IR luminosity density data at z = 0.0 and z = 0.1 (see 
i) 13. ip . These values can be found in Table [T] With the 
stellar component of the luminosity density calculated 
from eqn. (j4|), the dust component of the luminosity 
density can be determined. Both the stellar and dust 
luminosity densities are used to model the EBL energy 
density. 

2.3. Extragalactic Background Light Model 

The comoving EBL energy density at the current epoch 
can be obtained by integrating the luminosity density, 
eqns. ^ and ([TT|) . 



euEBL{e; z = 0) 

f 

Jo 



2 2 



dzi 



e'jie';z^) 
(1 + zi) 



dN 
dedV 
dt^ 



dzi 



(12) 



where e j(e; z) — g j«*i'"s(gj 2:) + ej *{e; z) and e' = (1 + 
zi)e. If star formation and star emission ended at some 
redshift z > 0, then the EBL energy density observed 
today would be 



euEBLie;z) 



dzi 



e'j{e';z,) 
(l + ^i) 



dt^ 
dzi 



(13) 



This is the energy in the comoving frame per unit co- 
moving volume as a function of e, which is the photon 
energy in the comoving frame. The energy and volume in 
the comoving frame are quantities in coordinates which 
expand with the universe, and hence are the quantities 
which we observe today. 

However, observers in a galaxy at a redshift z > 
would have viewed a universe that is (l + z)^ smaller than 
now, and photons observed now with energy e would be 
observed by them to have proper energy ep = (1 + z)e. 
The proper energy density they would have observed is 
given by 



e-p UEBL,p{e-p] z)^ {1 + z)4 e uebl{(-] z) 



(1 + ^r 



where 



dzi 



e'j{e';z,) 
(1 + zi) 

1 + Zi 



dzi 



(14) 



This is the radiation field with which a 7-ray photon 
emitted at redshift z interacts on its way to Earth. 



TABLE 1 

Dust Parameters. See text for details. 



Component 


n 


fn 


r„ [K] 


e„ [10-9] 


Warm Large Grains 


1 


0.60 


40 


7 


Hot Small Grains 


2 


0.05 


70 


12 


PAHs 


3 


0.35 


450 


76 



Energy density can be converted to intensity in units 
of, e.g., nW m~^ sr~^ by 



e/(e;z) = — €Uebl{^;z) 
47r 



(15) 



which is convenient for comparing to EBL intensity mea- 
surements. 

3. NUMERICAL RESULTS 

In this section, we compute the luminosity densities 
and EBL energy densities from starlight and dust for 
various combin ations of SFRs and IMFs. Although the 
iSalpeteii ()1955l ) IMF (^(m) oc m~^-^^) is still in common 
use, it now seems that it overpredi c ts the number of low 
mass stars. iBaldrv k Gl azebroold ()2003[ ) fit local lumi- 
nosity density data to find a preferred IMF, which is simi- 
lar to a Salpeter IMF above m = 0.5 {£,{m) oc m"^-^) but 
flatter below this mass {S.i'm) oc to~^-^). The Salpeter 
A IMF is also often used, which is ^(m) oc m~^'^ be- 
low m = 0.5 ancl £(m) oc m~^'^^ above m = 0.5. 
iHopkins fc B eacoml ()2006f ) compiled a large amount of 
SFR data and fit them using various IMFs. To parame- 
terize their SFR fits they used the function of lCole et al.l 
(|2001h as well as their own piec ewise power-law func- 
tion, and found the Salpeter A and lBaldrv fc GlazebroolJ 
(j2003) IMFs were favored, under the assumption of a uni- 
versal IMF. 

For our models we use the same designations as RDF09 
which are as follows (see RDF09 for details on these SFRs 
and IMFs): 



Model A: iCole et all (|200l[ ) SFR and Salpeter A 
IMF 



i Model B: [ Hopkins fc Beacoml (|2006l) SFR using 
iCole et all ( 20011 ) parameterization and Salpeter A 
IMF 



■ Model a 'Hopkins & Beacoml (poM ) SFR us- 
ing Cole et al. (2001) parameterization and 
IBaldrv fc Cxia^rookl lj2003h IMF. 



Model D: IHopkins fc Beacoml (p006l ) SFR and 
Salpeter A IMF 



• Model E: Hopkins & 
IBaldrv fc Glazebrook 



BeaTOmI (|200l ) SFR and 

(2oo|)1mf 



Measurements of SFRs are generally only sensitive to 
the formation of massive stars, so that an IMF must be 
chosen to extrapolate to lower masses and determine the 
global SFR. Thus, the SFR and IMF are not independent 
of one another and a combination of both are needed to 
calculate luminosity densities and EBL energy densities. 
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3.1. Luminosity Density 

The luminosity density at a particular wavelength is 
measured by adding up the luminosities of the galaxies 
at that wavelength in a particular volume, and form- 
ing a luminosity fun ction using the 1/Vmax method (e.g., 
iSandage et al.|[l979l ). This involves finding the luminosi- 
ties within a redshift slice, which gives the maximum co- 
moving volume {Vmax) within which the galaxies could 
be found. Since the surveys are limited to a minimum 
observable magnitude, some extrapolation to lower lumi- 
nosities must be done, which involves a fit to the lumi- 
nosity function. In the optic al and UV , the lu minosity 
function is typically fit with a iSchechted ()1976( ) function 
which is 



4>{L) 



(j)^ (L/L^)'^ exp(-L/LH.) 



(16) 



where the function has three fit parameters, i*, and 
a. This function can be integrated to give the luminosity 
density, so that 



dL L (j){L) = 0*L,r(a + 2) 



(17) 



where r(a;) is the Gamma function. The Schechter func- 
tion, written in terms of absolute magnitude, M (not to 
be confused with stellar mass in § [2|) , is 

4>{M) = O.41n(lO)0JO°-''(*^-*'^*)("+i) (18) 

xexp(-10-"-4(*^-*-^*)) , 

which can then be integrated to give the luminosity den- 
sity in terms of AB magnitudes per volume. 



jj^i = _ 2.5 log [0,r(a + 2)] + C 



(19) 



where C is the conversion from observed magnitudes to 
AB magnitudes. The error in the luminosity density is 



OM 



\dM, 



2 



dj_ 



da 



(20) 



/ 2.5 T{a + 2) a^, 
Vln(lO) 

/ 2.5 



Vln(lO) 



-00 (a + 2) CTc 



where 4^o{x) is the digamma function (not to be confused 
with tpiz), the SFR, used in ii l2.ip . This can be converted 
to power per volume (in, e.g., W Mpc~^) by 



ly X 10" 



and 



V\n{10) JAf -34-10 



-2.5 



10^^ 



(21) 



(22) 



where v is the central frequency of the bandpass. 

Luminosity functions in the IR can be quite different 
from the optical, and are often fit with broken power- 
laws, 

J L < 

m = { ' r;(i-/3 > (23) 

L > 



(e.g.. iBabbedge et al.ll2006l) instead of a Schechter func- 
tion. Integrating this, then, to get the luminosity density 
gives 



ejt J d.L (j){L) = (jj^L^ \j2~ 



1 



2-13 



(24) 



assuming /3 > 2. Aside from direct observations, the UV 
luminosity density at high redshift can be probed with 
measu rements of the photoio nization rate from the Lya 
forest (jFaucher-Giguere et al. 2008a, b). 

In terms of our theoretical formalism, the luminosity 
density is dependent on the IMF, SFR, and dust param- 
eters, and can be calculated from eqns. (U) and (ITT|) . The 
results of our models, plotted with measurements of the 
luminosity density at various redshifts, can be seen in 
Figs. [5] and [31 Although the differing measurements do 
not entirely agree, they provide a consistent picture for 
the evolution of the luminosity density, particularly in 
the UV/optical and at low redshift. 

It appears that our Model C best fits the luminos- 
ity density data. Models C and E fit the data about 
equally well at most wavelengths, but model C fits the 
UV data at low redshift (especially z = 0.1) better than 
does model E. These models agree with the data quite 
well at low z, and less-well at high z; however. Model C 
still provides a better fit to the data at high z than do 
the other models. Measurements at higher redshifts are 
less certain, and less likely to agree with each other. In- 
terestingly, the luminosity densities at high z from Lya 
measurements (jFaucher-Giguere et al.|[2008b[ ) are about 
a factor of two higher than me a sureme nts from direct 
imaging. iFaucher-Giguere et al.l ()2008bD speculate that 
this may be due to a faint population of galaxies at high 
redshift which have avoided detection in deep surveys 
so far. Other possible explanations for a discrepancy at 
high z include a changing IMF, and a changing absorp- 
tion from dust; we assume both of these are independent 
of re dshift in our models. T he 24 /im luminosity density 
from lBabbedge et al.l (|2006( ) increases rapidly at high z, 
and is considerably above our models and other data. 

The higher redshift luminosity density contributes lit- 
tle to the local EBL. Thus such discrepancies have little 
effect on the direct EBL observations (ii l3.2p and absorp- 
tion of TeV 7-rays from nearby blazars (§ H]), although 
it could have implications for absorption of 7-rays from 
higher z objects observed with Fermi or AGILE, how- 
ever, even here generally only for the UV/optical. 

3.2. Extragalactic Background Light 

Direct measurement of the EBL energy density is diffi- 
cult due to conta mination from fo r egrou n d zodiacal and 



Galactic light ("Bern stein et al.l 2002t_ 'Goriian et al 
2000; iDwek fc AreiidlJ Il998t fCambresv et al 
2001; iWright fc Reesd boOOt iLevenson et all l2007t 



iHauser et al.l I1998D . For example, the foreground 
subtraction of the optical measureme nts of the EBL 
with the Hubble Space Telescope by 'Bernstein et alj 
(I2002D have been vi g orously debated^ (Mattila 20O3 
IBernstein etHI [200l IBernsteinI [20071 ). The most 
accurate IR-to-UV EBL measureme nts are the measure- 
ments of the far-IR with COBE (jFixsen et al.l Il998f) 
and BLAST (jMarsden et al.ll2009 f). Galaxy counts ma 



be used to estimate the EBL ([Madau fc Pozzetti .200 
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Fig. 2. — Data and models of the comoving luminosity density at various redshifts. Data a re black symbo ls with error bars, as follows: filled 
circles; Blanton et al. (2003); filled squares: Cole et al. (2001); filled diamonds: Kochan ek et al.l l|20011 ): filled upward-pointing triangles: 

tidavari et al. (2005); filled downward-pointing trian gles: Tresse et al. (2007); filled leftward-pointing triangles: Sawicki & Thompsq3 
006); filled riglitward-pointing triangles: Dahl en et al.l (1200 7): star: Smith et al. (2008); empty upward-pointing triangles: Magnelli et al] 
009); emp ty diamon ds: Takeuchi et al. (2006); empty downward-point ing tr iangles: Babbedge et al. (2006); empty leftward-pointing 
triangles: iHuang et al.l ( 12007 ); empty rightward-pointing triangles: Le Floe h et al.i (l2005i ): circle with cross inside: iFlores et al.l (|1999|) . 
Curves are our models as follows: solid red: Model A; Dashed double-dotted green: Model B; short dashed blue: Model C; long dashed 
violet: Model D; dashed-dotted cyan: Model E. 



Fazio et al.ll200l iDole et alJ [20061 : [Metcalfe et aL|[200l 
Papovich et al.ll2004[ ) but the crowded fields make this 
difficult, and the unknown number of unresolved sources 
results in a lower limit. Lcvcnson & Wright (2008) have 
used a Monte Carlo Markov chain analysis to correct for 
crowded fields in 3.6 /im galaxy counts with the Spitzer 
Space Telescope, resulting in significantly stronger lower 
limits. Their new technique, when applied to other 
wavelengths, could drastically improve estimates of the 
EBL. Still, the lower limits from galaxy count s, even 
the improved one bv iLevenson fc Wrightl (|2008f ). are a 
factor of a few below the direct measurements. This 
indicates either that a significant fraction of the EBL 
originates from extremely faint galaxies not resolved in 
number counts, or there is significantly more foreground 
light than assumed. Franccschini ct al. (2008) consider 
the latter to be more likely, based on the fact that the 
galaxy counts are performed in fields with quite deep 
exposures. 

The EBL photon density, calculated with our mod- 
els, can be seen in Fig. along with direct measure- 
ments, lower limits from galaxy counts, and upper lim- 



its from TeV blazar observations. Models C and E are 
very close to the galaxy count lower limits in the optical 
and UV, and in fact are bel ow th e near-IR lower lim- 
its from iLevcnson & Wright ()2008D and the far-IR lim- 
its from [Dole et al.l (|2006). Models A, B, and D, on the 
other hand, seem to be marginally a bove the upper limits 
from 7-ray observations of blazars (|Mazin fc Rauel[2Q07l: 
iFinke fc Razzaau^ I2009() although these results a re not 
without controversv (e.s... iStecker fc Scullvll2008[ ). Due 
to how well Model C fit the luminosity density data, we 
take this as our best fit model. 

The integ rated intensity from our models can be seen 
in Table [H IFardal et al.l ()2007l ) created a simple model 
by a fit which compromises between direct photometric 
observations and number-count lower limits on the local 
EBL intensity. They find an integrated EBL intensity 
of 50 - 130 nW m~^ sr~^, and a value closer to the 
lower end is required to be consistent with the K-hand 
lum inosity density (s i milar to our best fit model, model 
C). INagamine et alJ (|2006D find this value to be 43 ± 
7 nW m~^ sr~^, which is in good agreement with our 
model C. 
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Fig. 3. — Data and models of the comoving luminosity density at various redshifts. Symbols and curves are the same as Fig. (2] while 
additional data are from the following sources: cross: [Cirasuolo et al. (2007); empty circles: Faucher-Gigucre et al. (2008a b); empty square: 
iRe ddv ct al. (2008); diamonds with horizontal stripes: Caputi et al. (2007). box with horizontal stripes: Perez-Gonzalez et al. (2005)- 



Our model of the steUar component is more accurate 
at low redshifts than high redshifts, and more accurate 
than the dust component, since there are more luminos- 
ity density data in this region with which to normalize 
our results. Results from our best fit model are in quite 
good agreement with lower limits from galaxy counts in 
the near-IR through UV, except for the recent 3.6 /xm 
point of Levenson & Wright (200l). In the far IR, this 
model is slightly below the points of iDole et a l. (2006), 
altho ugh they do fit the BLAST data (jMarsden et all 
120091 ) fairly well. The addition of a separate component 
representing ultraluminous infrared galaxies, as done by 
[Kneiske et al. (2004), may bring our model C in agree- 
ment with these other lower limits. However, it is un- 
clear why this component would not be represented in 
the luminosity density measurements (Figs. [2] and [3]). 
iKrennrich et al.l (12005 ) have recently used the lower limit 
of lLevenson h Wrighti ((200l) to construct their own EBL 
SED, which gives a greater energy density than our 
model C. 

Fig. [5] shows the redshift evolution of model C. This 
peaks at around z « 2, which is about where the SFR 
peaks (see RDF09). Also note that at higher redshifts, 
there are more high mass stars to create a larger high 
energy component. The high z SEDs are flatter in 
the optical-UV due to this relatively higher contribu- 



tion of high mass stars. At lower z, the contribution 
is greater from longer living low mass stars, as well as 
those high mass stars which have evolved off the main 
sequence. This also affects the dust emission compo- 
nent. At around z = \ the far-IR dust component peaks 
at a greater energy density than does the optical stellar 
component, and the mid-IR PAH component gets pro- 
gressively greater at higher z. Since the dust absorption 
is greater at higher energy (lower wavelength) , where the 
high mass stars emit most of their radiation, the absorp- 
tion is greater at higher z, and thus the dust emission is 
greater at higher z. 

The Balloon-borne infrared telescope BLAST found 
that the fraction of galaxies with redshift z < 1.2 which 
contribute to the local EBL is fsBLiz < 1.2, 250 fj,m) = 
0.60 ± 0.11, fEBLiz < 1.2, 350 ^im) = 0.49 ± 0.10, and 
fEBL{z < 1.2,500 fj,m) = 0.39 ± 0.11, where the error 
bars include measure ment and calibration uncertainties 
(jMarsden et al.l [2009D . We calculate these fractions for 
our model C, to test its consistency with these BLAST 
results, with 



/ebl(< z,e) = 



lo dzi 



e' J"(e';zi) 
(l+^i) 



(25) 



eUEBLii]Z = 0) 

Results from this calculation can be seen in Fig. |6l As 
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TABLE 2 

Integrated EBL intensity from our models, in nW 



Model 


Stellar Component 


Dust Component 


Total 


A 


36.0 


25.9 


61.9 


B 


36.4 


25.5 


61.8 


C 


26.5 


20.4 


46.8 


D 


37.0 


26.1 


63.1 


E 


25.5 


19.6 


45.0 



can be seen, these fractions are in excellent agreement 
with the results from BLAST. 

3.3. Comparison with other models 

We have plotted our best fit model, Model C, along 
with several other models from the literature in Fig. [71 
N ote that our model agree s re markably well with th ose 
of lFranceschini et all (|2008D and lGilmore et al.l ()2009[ ) at 
> 0.3 eV (< 5 /zm). Of t he models pres e nted in Fig. [3 
only the best fit model of Kneiske et al. (2004) and the 
fast evolution model of Stacker e t ah (.20 06) agr ee with 
the re cent 3.6 /im lower limit from iLevenson fc Wright] 
((20081) . 

Two previous models from RDF09 are plotted here, us- 
ing the single power-law version of the stellar luminosity- 
mass relation. These previous versions of our model in- 
clude the contribution to the EBL from main-sequence 
stars. Model B from RDF09 was considered our best fit 
model from this previous work, while we have plotted 
the old Model C for comparison with our best fit model 
from this work. It is clear that post-main sequence stars 
have a considerable effect on the result, as the C models 
differ by as much as an order of magntitude at longer 
wavelengths (■^ 1 — 10 fim). The model in this work is 
more intense longward of '~ 1 fim due mainly to giant 
stars, while shortward of ~ 0.3 fim the inclusion of white 
dwarfs explains the discrepancy between these works. 

The modeling technique of Kneiske et al. (2002, 200|) 
is quite similar to ours: they also integrate over stel- 
lar parameters and treat dust extinction and re-emission 
in a similar fashion. Our results, however, differ sig- 
nificantly from theirs due to our use of updated pa- 
rameters for the SFR, IMF and dus t extinction. The 
best fit model of 'Knei ske et al.l (j2004) also includes sep- 
arate components from ultraluminous infrared galaxies 
(ULIGs) and normal galaxies, which, if applied to our 
model (as mentioned above) , could result in a local EBL 
energ y densities i nore in ag reement with the lower lirn .- 
its of (Dole et a T (2006') andE evenson fc Wrightl (l200l . 
with which the Kncisk e"et al.l ()2004D best fit model is 
in agreement. Their infrared component, however, con- 
siderable ovBr2redictstlie_^0 /im measur ements from 
BLA ST (|Marsden et al.l [20091) and COBE ([Fixsen et alj 
|1998() . Disentangling the dust extinctio n from these 
two di fferent p opulations, using the rea listic [Driver et alJ 
^0§) (or the iPopescu fc filffl [2009 ^ model would be 
considerably more complicated. Furthermore, greater 
dust emission may contradict the luminosity density data 
at low z (Fig. [1 

iStecker et alT(|2006D have used a backwards evolution 
model based on observations of the infrared luminos- 
ity function to calculate the IR EBL energy density. 



They then use a f orwar d evolution model similar to 
iSalamon fc Steckeil ([19981 ) to calculate the optical/UV 
EBL, and normalize it to their backward-evolution cal- 
culation. Our model peaks in the IR at an intensity 
relatively close to their baseli ne model, both o f which 
are below the IR lower limits of [Dole et al.l ([2006[ ) . How- 
ever, our models differ by as much as a factor of 8 in 
the UV. This could be because extrapolating below 10 
/im may not be appropriate for de termining the infrared 
background ([Lagache et al.l l2003f ) . This wavelength is 
critically important to the models of[S"tccker et al. (2006|) 
because this is where they normalize the results of their 
forward-evolution modeling, and the optical/UV region 
is the most important for 7-ray absorption (§[1]). Further- 
more, their model has neglected dust extinction, which 
has a considerable effect at UV wavelengths. We note 
that the far IR peak of their baseline model, which is 
based in IR luminosity density data, is within 10% of our 
far-IR peak, which is fit to similar IR luminosity density 
data. 

iGilmore et al.l ([2009[ ) have presented models based on 
semi-analytic models of galaxy evolution, using cosmo- 
logical parameters from WMAP and normalizing their 
semi-analytic model's luminosity density prediction to lo- 
cal luminosity density measurements. They also include 
a component from quasars and realistically treat extinc- 
tion of ionizing radiation by the intergalactic medium, 
both of which have implicat ions for the UV port ion of 
the EBL at high redshift. IGilmore et "aD ((2009) only 
treat the UV-optical portion of the EBL, while their 
further work will be presented in Gilmore et al. (in 
preparat ion). This is an update o f prev ious, similar 
models ([Primack et al.l [19991 [200l [20051) . The stel- 
lar component of our model probably agrees with the 
IGilmore et al.l (2009) fiducial model so well because we 
have both made sure our models agree well with luminos- 
ity density measurements (§[33J- However, at lower ener- 
gies (longer wavelengths) their model intensity is as much 
as double ours. Moreover, their predicted inten sity is 
above the 250 fim BLAST measurement (Marsden et al.l 
'2009") , although it is consistent with the COBE results 
((Fixsen et al.lll998[). 

The model of iFranceschini et all ((2008() is based on 
integrating luminosity functions at various wavelengths 
and redshifts to generate an EBL energy density. Again, 
their model agrees well with ours at higher energies 
but at 20 /im they predict about three times the in- 
tensity as we do. The infrared component of their 
model is comparabl e to that of the best fit model of 
iKneiske et al.l (l2004f). and similarly, i t over-predicts the 
250 /im BLAST fMarsden e t iI[|2009D measurement and 
the COBE (Fixsen et al. 199^infrared data. 

4. ABSORPTION OF 7-RAYS 

Once the energy density of the EBL has been calcu- 
lated using eqn. ([Ti| . the absorption optical depth of 
7-ray photons as a function of observed 7-ray photon 
energy, ei, can be calculated by 



y(ei,z) 



cnr^ 



dz' 



dtf, 



dz' 



ipUEBL,p 



(1 + z')' 
0(epei(l + z')) , (26) 



1(1+^') 
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Fig. 4. — Our models for the EBL plotted along with measurements and cons traints from o bservations. The cur ves are our models 
with the same symbols as Fig. |2] Measurements arc from 'Bernste in et al.l l|2002!, cyan points), 'Gorjian et al.! f200(J, empty red circle^, 
[D^ck & Arcndt (1998, green asterisk), Cambrcsy et al. (2001, empt y cyan squ are), Wright & Reese (2000, black cross), Lcvcnson ct aU 
^2007, maroon diamonds), Hauscr et al. (1998, green filled circles), Fixsen et al. (1998, brown filled diamonds), and Marsden et al. (2003^ 
violet filled squares). Lower limits arc from Fazio et al. (2004, red empty triangles), Madau & Pozzetti (2(500, brown filled triangles), 
[Eevenson & Wright (2008, blue filled triangle). .Dole ct al.. (2006; . ma genta filled triangles) , .Metcalfe ct al. (200^, black em pty triangle), and 
[Papovich ct al. (2004, green empty triangle). Upper limits are from iHauser et aLl ( 119981 . brown filled inverted triangles). IDwek fc Arendtl 



(ir998i . blue empty inverted triangles), IMazin &: Hauel ( 120071 . upper and lower black curves F^*^" = 0.67 and - 
respectively), and red empty and black filled inverted triangles are the F™^" = 1.0 and F™j" = 1.5 upper limits, 
ll^'inke &: Hazzaquel (120091 ). The black curve at long wavelengths is the cosmic microwave background. 
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Fig. 5. — The proper EBL energy density as a function of proper 
photon energy for model C, for a variety of redshifts. 




Fig. 6. — The fraction of the local (2 = 0) EBL which originates 
below a redshift z at the three BLAST wavebands ijMarsden et al.l 
[2009D for our model C. 



where \dt^/ dz\ is given by eqn. ([S]), is the proper frame 
EBL photon energy, uF^g]^{e^] z') is the proper frame 
EBL energy density given by eqn. (|14p . 

-^(so) = T^^lnwo -/3olnwo - -j-^^ (27) 

1 Pg -L PO 

+2/3o + 41nwoln(l + M;o) -4L(wo) , 



[il = l- 1/so, wo = {l + /3o)/(l - /3o), and 

/too 
dww-^\n{l + w) (28) 

(|Gould fc Schredei<[l967l: IBrown et aLlfTOTl . 

We have used our m odel C to deabsor b seve ral of 
the blazars presented in iFinke fc Razzaaui ()2009D . We 
choose blazars that have a low ^ parameter, defined in 
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Fig. 7. — Our model C (solid black curve), the best fit 
mo del of IKneiske et alj II 2004I . red dotted curve) , the model of 
IFranccschini ct al. (2008, short dashed green curve), the fiducial 
model of Gilmorc ot al. (2009, long das hed blue curve) and the 
fast evolution and baseline model from Stec ker et al.l l|200g . up- 
per and lower dot-dashed violet curves, respectively). The double 
dot-dashed and the dot double-dashed black curves are the single 
power-law model B and model C from RDF09, respectively. 

iFinke fc Razzaaud (|2009D to be 

^ ^ f ^"'^ ~ " In f^IIl^\ ^ (29) 

(not to be confused with ^(m), the IMF in § 12. 1|) where 
Erain and Emax are the minimum and maximum energies 
of the VHE 7-ray spectrum, Fobs is the observed photon 
spectral index, and F™™ is an assumed minimum intrin- 
sic photon spectral index, for which we use Fi„t — 1.5, 
the lowest value one would expect from naive test parti- 
cle shock acceleration theory where particles are acceler- 
ate d with number indices stee per than —2 (for a review, 
seelBlandford fc EichleiifTOSTl) although see lStecker et al.l 
(|2007ir for a different interpretation. The lower the pa- 
rameter ^, the lower the deabsorbed spectral index F 
should be. We also choose the VHE 7-ray spectra of two 
bright, distant quasars, 3C 66A and 3C 279 to deabsorb, 
to give a greater range in redshift. The deabsorbed spec- 
tra are fit with power-laws, and the results can be seen in 
Fig. [5] and Tabled Note that in Fig. [5] the vF^, spectra is 
plotted, where F = 1.5 corresponds to vF^ oc E'^-'° . The 
results are generally consistent with a minimum intrinsic 
Vint of 1.5, as this value is within the error bars for all 
of the sources. 

A plot of the energy at which the universe becomes 
optically thick to 7-rays, defined where r^^ = 1 is given 
in Fi g.[9l which is often known as the Fazio-Stecker rela- 
tion (jFazio fc~S teckcr 1970), for several models including 
our Model C. Also plotted are the maximum photon en- 
ergy bin of several blaza rs, observed with atmospheric 
Cherenkov telescopes fsee lFinke fc R azzaque 200 9j, for a 
list aii d references) an d the GRB s 080916C (,Abdo et al.l 
I2009b[) and 090902B ()Abdo et a l. 2009a) observed with 
the i^ermi-LAT. The VHE 7-rays from many blazars are 
highly attenuated by the EBL, since several are consid- 
erably above the r^^ = 1 for all models, and the highest 
energy photons from GRBs constrain the EBL at high 
redshifts. Also note that the universe will be optically 
thin to 20 GeV and lower pho tons over all red shifts for 
all models except those of iStecker et al.l (|2006() . 



We have plotted the absorption optical depth for our 
model and several others in Fig. [TOl The models are 
optically thin (r^-y < 1) at 200 G eV fo r z < 3, ex- 
cept the models of IStecker et al.l (j2006[ ). Observing 
a large number of high-energy photons with observa- 
tories such as Fermi and AGILE from high-z sources 
could rule out these models. Our model is moderately 
more opaque than those of iGilmore et al.l (|2Q09f ) and 
iFranceschini et all ()2008D at these redshifts, although it 
would be difficult to distinguish these models based on 
7-ray observations. 

5. SUMMARY AND CONCLUSIONS 

We have developed a model for the UV through far-IR 
EBL from direct stellar radiation and stellar radiation 
which is absorbed and reradiated by dust. This approach 
extends RDF09, which applied between ~ 1-10 eV. Our 
best fit model. Model C, is consistent with the collec- 
tion of measures of SFR bv iHopkins fc BeacomI (j2006), 
and a variety of luminosity density data at various red- 
shifts z < 3. It does not require complex stellar structure 
codes or semi-analytic models of galaxy formation. Our 
model is most accurate at low- 2; and high e (short wave- 
lengths), which is the energy and redshift range that is 
useful for calculating the absorption of TeV 7-rays from 
nearby blazars by electron-positron pair production. The 
EBL energy density of Model C is generally consistent 
with lower-limits from galaxy counts and is below direct 
EBL measurement s except for the lower li mits at 3.6 
/im (Lc venson fc Wright! [20081 ) and 60 /xm pole et al.l 
2006) . It also agrees w ell with the recent EBL mod els 
of .Gilmore et al.l (|2008D and IFranceschini et al.l (|2008[ ) in 
the near-IR through UV, although less well at longer 
wavelengt hs. Our Model C is significantly below the 
models of IKneiske e t al. (2004) and Stecker et al. (200i) 
for all wavelengths. 

The SFR, IMF, and dust parameters in our best fit 
EBL model were chosen to agree well with luminosity 
density measurements, particularly those at low-z where 
the data is best, and to give precedence to luminos- 
ity density measurements over local EBL measurements 
when deciding the quality of an EBL model. Although 
different luminosity density measurements are not always 
consistent with each other, as one can see in Figs. [Hand 
[3l they are less controversial than m easurements of the 
local EBL intensity. For example, IMadau fc Pozzet"t]| 
((2000.) claim that their galaxy count m easurements re- 
solve almost all of the flux in the EBL. [Bernstein et al.l 
(|2002D suggest, however, that those measurements fail to 
take into account the faint portions of t he distant galax- 
ies. T he more direct measurements of [Bernstein et ahl 
(|2002l ) give higher values than galaxy counts, but their 
measurements have also been criticized for the sub- 
traction of zodiacal foreground light (e.g. IMattilal 120031 : 
[Bernstein et al.l2005HBernsteinll2q07D. To reconcile these 
differing inferences, " iTotani et al.l (|200H) suggest the pos- 
sibility that this discrepancy can be resolved by a faint 
EBL component outside of normal galaxy populations. 

By constructing an EBL model consistent with 
luminosity-density observations, we have found that 
an EBL energy density very close to the lower lim- 
its from number cou nts is require d. Similar conclu- 
sions were reached bv i Fardal et al.l (|2007D . They point 
out that their results are strongly dependent on the 
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Fig. 8. — TeV blazar spectra observed (circles) and deabsorbed by our model C (squares). 



TABLE 3 
Deabsorbed TeV Blazars 



Blazar 


Redshift 


Observed F 


Deabsorbed F 




lES 2344+514 


0.044 


2.54 


± 


0.18 


2.10 


± 


0.68 


4.7 


lES 1959+650 


0.047 


3.18 


± 


0.17 


2.18 


± 


0.33 


4.7 


lES 1426+428 


0.129 


2.60 


± 


0.60 


1.91 


± 


0.58 


2.2 


lES 0229+200 


0.139 


3.09 


± 


0.26 


1.07 


± 


0.45 


1.8 


lES 1101-232 


0.186 


2.94 


± 


0.20 


1.27 


± 


0.46 


5.7 


lES 0347-121 


0.188 


3.10 


± 


0.25 


1.56 


± 


0.43 


6.5 


3C 66A 


0.444=' 


4.10 


± 


0.72 


1.28 


± 


0.98 


19.6 


3C 279 


0.536" 


4.11 


± 


0.68 


2.33 


± 


0.89 


16.0 



Reference 



Schrocdter et al 
Aharonian et al. 
Aharonian et al. 
Aharonian et al. 
Aharonian et al. 
Aharonian et al. 
Acciari et al. (2009 1 
Albert et aT"(2058l 



(2005) 
(2003a) 
(2003b') 
(2007c) 
(20q7a) 
(2007b) 



IMF, which may vary over the history of the universe. 
Although radiation hydrodynamic simulations indicate 
that radiative feedbac k could lea d to a universal IMF 
in the local universe (jBatd [20091 ) . a detailed modeling 
of a large amount of Sloan Digital Sky Survey data 
indicates that fainter galaxies generally have an IMF 
whi ch produces fewer massive star s than brighter galax- 
ies ([Hoversten &: Glazebrookir2008[ ). This result has fur- 
ther evidence in t he ratio of Ha to f ar- UV flux in H I- 
selected galaxies (jMeurer et al.ll2009| l. iBabbedge et alJ 
(j2006) suggested that a non-universal IMF may be re- 
sponsible for their measured rapid increase in the 24 fim 

l uminosity density fi) 13.11) . 

iGeorganopoulos et al.l (|2008D have suggested a new 



method for measuring the EBL. They point out that 
Compton scattering of the infrared component of the 
EBL by the radio lobes of Fornax A could be detectable 
by Fermi in ~ 2 years in scanning mode, and that it 
may also be possible to put upper limits on the optical 
EBL component in the absence of Fermi detection. This 
will be particularly important for distinguishing our best 
fit model from those which pr edict significan t ly hig her 
IR emission, such as th ose of iKneiske et al.l (|2004l ) or 
iFranceschini et al.l (|2008l ). 

In the near future, surveys such as the Dark En- 
ergy Survey, and the Sloan Digital Sky Survey II will 
detect ~ 10'^ core-collapse supernovae per year, and 
Pan-STARRS and the Large Synoptic Survey Telescope 
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Fig. 9. — A plot of T-j-y = 1 for several EBL models including our 
model C, as a function of redsliift. Curves signify the same models 
as in Fig. [7] Also plotted by the filled diamonds are the maximum 
7-ray photon energy bins from blazars observed with atmospheric 
Cherenkov telescopes and the GRBs 080916C and 090902B ob- 
served with the Fermi-LAT. 

will detect ""^ 5 x 10 ^ core-collapse supernovae per year 
(|Lien fc Fields! |2009[) . This will lead to extremely pre- 
cise measurement of the SFR, as the core collapse su- 
pernovae rate is a good tracer of the high mass SFR, 
although modeling dust attenuation will still be a sig- 
nificant challenge. This in turn will lead to strong con- 
straints on models of the EBL and luminosity density. 



Detection of the neutrino background from core-collapse 
supernovae wil l furth er constrain the SFR and EBL (e.g., 
iHoriuchi et al.l 120081 ). Combinations of these various 
measurements mean that very soon the local EBL will 
be known to much greater precision. 

We have used our best fit EBL model (Model C) to cal- 
culate the absorption optical depth to 7-rays from cosmo- 
logical sources, valid over all relevant 7-ray energy ranges 
and redshifts. Our results show that absorption is quite 
significant for blazars observed at TeV energies by at- 
mospheric Cherenkov telescopes, and that de-absorbed 
VHE blazar spectra give results generally in agreement 
with particle acceleration theory. We also find that the 
universe is transparent (x^^ < 1) for all redshifts at en- 
ergies less than 20 GeV, which are those most relevant 
to the Fermi-LATQ. 



We are grateful to the anonymous referee for useful 
comments which have improved this work, including the 
discovery of an error in the original version of Fig. 5. 
We thank P. Eggleton and C. Tout for correspondence 
regarding corrections to their stellar formulae, and A. 
Franceschini, R. Gilmore, and T. Kneiske for correspon- 
dence related to their EBL models. This work was sup- 
ported by the Office of Naval Research and GLAST Sci- 
ence Investigation DPR-S-1563-Y, and by NASA Swift 
Guest Investigator Grant DPR-NNG05EE)41I. 
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